Coherence Properties of the Microcavity Polariton Condensate 
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A theoretical model is presented which explains the dominant decoherence process in a microcav- 
ity polariton condensate. The mechanism which is invoked is the effect of self-phase modulation, 
whereby interactions transform polariton number fluctuations into random energy variations. The 
model shows that the phase coherence decay, g™(T), has a Kubo form, which can be Gaussian or 
exponential, depending on whether the number fluctuations are slow or fast. This fluctuation rate 
also determines the decay time of the intensity correlation function, g^ 2 '(r), so it can be directly 
determined experimentally. The model explains recent experimental measurements of a relatively 
fast Gaussian decay for g^\r), but also predicts a regime, further above threshold, where the decay 
is much slower. 

PACS numbers: 71.36.+C, 78.20.Bh, 42.50.Lc 
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Microcavity polaritons arc quasi-particles arising from 
the strong coupling between excitons and photons con- 
fined in planar cavity structures. The observation of co- 
herent emission from a CdTe microcavity [l[ has demon- 
strated that polaritons can form a new type of quantum 
condensate. As in other quantum condensates, such as 
atomic gases or superconductors, a key property is the 
existence of an order parameter, the local phase, which 
is correlated over large times and distances. The polari- 
ton condensate presents interesting theoretical challenges 
since, unlike these other systems, it is mesoscopic, typ- 
ically consisting of a few hundred particles, and out of 
equilibrium, with pumping required to maintain the pop- 
ulation against emission losses. In mesoscopic systems, 
order parameters fluctuate 0, Q , so the phase correla- 
tions decay. In this letter, we present a theory which 
shows that the source of the fluctuations is variations in 
the number of condensed particles, induced by the pump- 
ing process: it is this dynamics that is responsible for 
the decoherence. Our theory shows that, under appro- 
priate pumping conditions, existing microcavity struc- 
tures should display much longer coherence times than 
currently measured, opening up opportunities for exper- 
iments manipulating the quantum state of the system. 

The coherence of the polariton condensate can be 
quantified by the decay of the first-order coherence func- 
tion, g^'ir) oc (a^(O)a(r)), whose Fourier transform is 
the emission spectrum. For the polariton condensate, 
this function is directly revealed by coherence measure- 
ments on the optical emission [H 0, 0| ■ In a condensate we 
expect long phase-coherence times and thus a spectrally 
narrow emission above threshold. Recent experiments 
[3| have shown that the decay time of ff^jMis ~ 150ps, 
much longer than was originally believed [J Q , but short 
compared to a laser or atomic gas. Furthermore, the de- 
cay has a distinctive Gaussian form. The experiments 
also determine the intensity-intensity correlation func- 
tion (second-order coherence function), <t 2 '(t), which re- 



veals significant number fluctuations (g^ 2 '(0) > I), de- 
caying with a timcscalc ~ lOOps. 

In this paper, we argue that the observed coherence 
properties of the polariton condensate are produced by 
the combination of polariton-polariton interactions, the 
number fluctuations generated by the pumping, and a 
critical slowing-down of these number fluctuations which 
occurs near the condensation threshold. One or more of 
these features are missingfrornprevious theories of con- 
densate coherence 0,0, i, Eft EH EH, yet all appear 
to be needed to obtain the behaviour now seen experi- 
mentally. Rather than attempting a detailed microscopic 
model, our treatment is based on general considerations 
of interacting, open condensates. Since our predictions 
agree well with the experiment, the model may be a 
useful guide in the development of microscopic theories 
@, @, [1, 0, [ll, El- The agreement with the experi- 
ments was briefly described with the experimental results 
in Ref. [1], focusing on the case of slow number fluctua- 
tions appropriate for the available data. Here we provide 
a comprehensive treatment of the theory, and show it 
predicts that other types of behaviour may also be ac- 
cessible in the experiments, most notably much longer 
coherence times due to motional narrowing. 

We first explain how the Kubo stochastic line-shape 
theory [3, [lj| can be applied to give a general form 
for g^'(r), accounting for both interactions and number 
fluctuations. This expression reduces to a Gaussian decay 
in the limit where number fluctuations decay slowly. We 
then develop a solution to a quantum-mechanical model 
of a pumped condensate, which gives the same the Kubo 
form for g^(r). Close to threshold, this solution gives a 
slow decay for number fluctuations, and hence explains 
the observed Gaussian lineshape, as well as the observed 
decay time. However, it also predicts that further above 
threshold number fluctuations will be much faster, and 
significantly longer coherence times should be obtained. 

In our discussions, we neglect spatial effects, which is 
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justified in the CdTe system where the emission spot 
is strongly localised by disorder. We thus model the 
condensate mode as a single anharmonic oscillator, with 
Hamiltonian 



H = a f a 



UJQ 



k (a 1 aY 



(1) 



where ujq is the oscillator frequency and k the strength 
of the polariton-polariton interactions. 

For a condensate of interacting particles, the inter- 
actions translate number fluctuations in the condensate 
into random changes to its energy, and so the coherence is 
lost. A similar effect, commonly termed 'self phase mod- 
ulation', was originally observed for laser beams propa- 
gating in a non- linear Kerr medium [l5l |. If we assume 
that the condensate has a Gaussian probability distribu- 
tion for the number of polaritons with variance cr 2 , it is 
straightforward to obtain g^ (Y) when the number fluc- 
tuations are sufficiently slow [4j . It has a Gaussian form 



\g^(r)\ = exp(-2/cW) = exp(-rV^) (2) 

As detailed in Ref.Q, we can obtain, directly from the 
measured data, values of er 2 ~ 25000 and k ~ 2 x 
10 _5 ps _1 , giving a decay time r c ~ 200ps, in reason- 
able agreement with the experiments. The value for k is 
consistent with theoretical estimates ljl 17 1 of the inter- 
action in a mode of linear size ~ 

The picture described above is essentially static; it as- 
sumes that the time-scale on which the number of polari- 
tons changes, r r , is much longer than the coherence time 
r c , so the only relevant time evolution is caused by the 
action of the Hamiltonian. The obvious problem with 
this description is that the coherence time, r c , is much 
longer than polariton life-time, tq ~ lps, due to emis- 
sion from the cavity. This suggests that the microcavity 
system may well not be in the quasi-static regime, and 
we need to consider the processes by which the number 
fluctuations occur in more detail. 

The effect of introducing a time-scale, r r , for fluctua- 
tions can be understood using the Kubo stochastic line- 
shape theory[3], which describes the decay of g^(r) for 
emission associated with a transition whose energy varies 
randomly in time. In our case, the random variations 
in energy are a consequence of the number fluctuations, 
with time scale, r r , determined by the pumping. When 
the random energies have a Gaussian distribution, the 
Kubo theory gives 



5 (1) (t)| = exp 
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(3) 



clearly be in the motional narrowing regime, giving an 
exponential function, with a very slow decay ~ 10ns. The 
measurement of the ff^fV) decay thus shows that r r is, in 
fact, much longer than To, and must be comparable to, or 
greater than, r c . This prediction is very well confirmed by 
the measurement of the decay of (r) , which directly 
determines the time scale for number fluctuations; the 
experimental decay time is ~ lOOps, similar to t c for the 
same condensate population. 

The explanation for the slow decay of </ 2 '( T ) comes 
from laser physics, where it is well known that number 
fluctuations are slowed close to the threshold. This can 
be explained using a simple classical model where the 
pumping provides a gain which is saturable, that is, has 
a dependence on the mode population. Above threshold, 
the linear part of the gain term exceeds the loss, so the 
population grows, and a non- linear saturation term is 
required to obtain a finite steady state. However, the 
response to small fluctuations, which is what the intensity 
correlation experiment measures, depends only on the net 
linear gain. Thus, near threshold, where the linear gain 
and loss are closely matched, fluctuations in the system 
relax very slowly. Haken[l|| shows the close analogy of 
this behaviour to the critical slowing down of fluctuations 
in the vicinity of an equilibrium phase transition. 

To put these considerations on a more formal footing, 
we have developed a quantum model of the polariton con- 
densate which can be solved analytically for g^(r) and 
g^ir). This model is a generalisation of one studied 
by Thomsen and Wiseman [7j, in the context of atom 
lasers, to include the regime where motional narrowing 
occurs; their model only applies to the 'far above thresh- 
old' limit, where the slowing of number fluctuations no 
longer occurs. The coherent mode is treated as an anhar- 
monic oscillator with a Kerr non-linearity, Eq. {T]). This 
mode is coupled to a reservoir, using the master equation 
formalism for the density matrix p. Reservoir losses are 
offset by a standard laser-like saturable pump term [l9| . 
We thus obtain equations for the population distribution, 
Pn = Pn,n, and the coherence u„ = \/np n -i,n 
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and 
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For T r slow compared to the scale set by the variance 
of the random energy distribution (r c above) , we are in 
the static regime and a Gaussian decay with life time 
r c is predicted. However, in the opposite regime, r r <C 
r c , motional narrowing occurs, and the decay becomes 
a simple exponential, with life time r ( 2 /2r r , much longer 
than r c . If we naively take r r = tq ~ 2ps, we would 



Here, ^7 = 1/tq is the cavity decay rate, and n c and n s 
are parameters describing the pump process plj. n c char- 
acterizes the strength of the pumping, while n s provides 
the saturation: for n > n s the gain decreases. Physically, 
this corresponds to the depletion of the pump reservoir 
by the processes which populate the condensate. Far 
above threshold, where the mean occupation (n) n s , 
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FIG. 1: (a) The decay of g^'(r) for populations (n) =50, 
500 and 1000. On this logarithmic plot, a simple exponential 
is linear, while a Gaussian is quadratic, (b) Decay times for 
g w (r) and # (2) (Y), as a function of population, obtained from 
the numerical solutions of Eqs.([SJ and @ (solid lines). The 
nonlinearity is k = 4 x 10 _5 7, and n s = 2.5 x 10 4 , as derived 
from the experimental results. The dashed lines are from 
the analytic expressions, Eqs.© and @. The marked points 
indicate the values of (n) used in (a). 



this model becomes identical to Ref. Q . Pumping terms 
like this, at least to 0(n 2 ), are required to give a finite 
condensate population, so the forms (0J and ([5]), with 
quantitative n c and n s , should be derivable from micro- 
scopic treatments of the polariton relaxation dynamics 
using kinetic equations [101 ] . 

The mean field threshold for condensation corresponds 
to n c = n s ; for n c > n s , the steady state solution of 
Eq.® is 



(n + n s )\ 



exp 



{n-{n)f 



2n c 



(0) 



where the Gaussian form is valid when n = n c — n s 3> 
1. The variance of the population is a 2 = n c , so that 
number fluctuations in the threshold region, n <C n s , are 
super-Poissonian. It is convenient to divide the threshold 
region into two: for n > 3yfnl ~ 500, the Gaussian form 
is valid for all n, because the non-physical n < states 
are not significantly occupied, while for smaller values 
of n, these states have to be explicitly excluded. In the 
former case, the mean population (n) = n. 

Eqs.® and ([5]) are easily solved numerically for pop- 
ulations of a few hundred particles. Fig. J!})) shows the 
decay times for the correlation functions obtained from 
these solutions, plotted as a function of the mean popu- 
lation. As the population increases, the g^ decay time 
rises rapidly to ~ lOOps, then decreases. The g^ 1 ' time 
also rises quickly, then flattens for a while before increas- 
ing again. Populations up to ~ 500 correspond to the 



experimental regime of Ref.0], and the observed g^ be- 
haviour has a very similar form. In Fig.flJ,), the actual 
decay of g^(r) is plotted, for three population values. 
For (n) = 50, corresponding to the initial rise in co- 
herence time, the form is a simple exponential. In the 
flat region, with (n) = 500 the decay starts off Gaus- 
sian, before becoming exponential at longer delays; this 
is the near static behaviour of the self-phase modulation 
regime. The final rise occurs when the g^ time short- 
ens and motional narrowing sets in. This is evident in 
the decay curve for (n) = 1000, which starts off with the 
same Gaussian as (n) = 500, but much sooner becomes 
a slower exponential. 

We now turn to deriving analytic solutions to Eqs.® 
and ([5]), which provide the dashed lines on the fig- 
ure. These solutions are valid in the regime where 
n > 3y/n^ ~ 500 and the Gaussian in Eq.® is valid 
without truncating the n < part. The solutions in- 
clude a term of include a term of exactly the Kubo form 
Eq.®, showing that this is the correct description of the 
physics. 

Eq. ((3]) for the population is independent of the non- 
linearity, so we can quote standard laser theory re- 
sults [191 ] for the intensity correlation function: 



9 {2} (r) = l + 



n s 

=2 eX P 



n \ n s — 

— 7T = l + =2 cxp(-7r), (7) 
n, J 7T 



where 7 = rvy/n c is the slowed decay rate. This result fits 
the single experimental data point fairly well; using the 
experimental measurements g( 2 ' (0) = 1.1 and n — 500 we 
obtain n s = 25000 and 7 = 7/5O. This gives a decay time 
of ~ 50ps, in reasonable agreement with the measured 
lOOps. 

The first order correlation function, 3^(t) cx 
J2 n u n(T)e luJ ° T , is obtained by solving Eq. The so- 
lution is required with an initial condition u n (0) = nP^ , 
which is a similar Gaussian function to but with 
mean n + n c /n. To obtain the correlation function, wc 
follow the approach of Gardiner and Zoller Using 
a Kramcrs-Moyal expansion, the difference operators are 
converted into differentials, leading to a Fokkcr-Planck 
equation for u, which we now write as it(n, t), with n a 
continuous variable. To deal with the appearance of n 
in the denominator of the pumping terms, we linearise 
around the mean value, writing n = n + n c /n + n'. Thus 
we obtain 
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where we have omitted constant nonlinear contributions 
to the oscillator energy, which can be absorbed in a 
renormalised luq. This equation is solved in the Fourier 
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domain [20( to give 



gW(r)\ = exp 



(e 7 +7T- 1) 

7 



x cxp 



(9) 



The first factor is just the Kubo expression, Eq. ([3]), with 
t 2 = 1/2k 2 <j 2 and <r 2 = n c , as before, and r r = 1/7. This 
constitutes the main result of our treatment: we obtain 
Kubo type behaviour, with the fluctuation time t t given 
by the decay time of g^ 2 \r). The second term corre- 
sponds to the Schawlow-Townes decay, enhanced in the 
threshold region due to the finite amplitude fluctuations. 
It is generally much slower than the first, in the regime 
where the expression is valid. 

The figure shows very clearly the importance of the 
fluctuation time scale r r oc 1/n (for large n) on the co- 
herence time. When n is increased, r c only changes very 
slightly, in fact decreasing as n c = n s +n grows, but r r 
shortens rapidly. This pushes the system into the mo- 
tional narrowing regime, where the life time T 2 /2r r is 
proportional to n. Note that we have to be careful treat- 
ing this as a prediction of a linear relationship between 
coherence time and emission intensity for high powers; 
we have simply kept n s constant and increased n c . In 
the textbook laser model n s is indeed independent 
of pump power P, and n c oc P. Though pumping of the 
the polariton system is considerably more complex than 
this, it is likely that n c increases more rapidly with P 
than n s does. Thus the fluctuation decay should decrease 
with stronger pumping, potentially taking the conden- 
sate into the motional narrowing regime. For this to be 
observable, it would, of course, be necessary that other 
mechanisms should not take over and restrict the coher- 
ence as the number fluctuation effect is suppressed. 

One surprising feature of the experiments is that very 
significant slowing down is still occurring at pump pow- 
ers of twice the threshold value, Ptn- This is inconsis- 
tent with the simplest assumption, n c oc P, since then 
7/7 = (1 — P-th/P) is close to one at these powers. A full 
understanding of this requires microscopic models of the 
pumping, but it may be explained by assuming that there 
are a limited number of reservoir states which can provide 
gain for the condensate, so the gain parameter n c must 
become independent of P at high pumping. In a system 
where the maximum achievable gain only just exceeds 
the loss, the reservoir must be almost full at threshold, 
so this saturation will extend the critical region to large 
values of P/Pth- 

Instead of increasing the pump power, an alternative 
way of obtaining larger populations and hence longer co- 
herence times may be to increase the size of the con- 
densate. In the experiments, the size is determined by 
the disorder, so this means finding larger emission spots. 
We present a simple argument for how the various time 
scales should vary with spot size, A, assuming constant 
density, that is n cx A. Since n = n c — n s , this suggests 



n c ~ A 1 and n s < A 1 , for large A. The polariton life- 
time, To is independent of A, while k ~ A^ 1 . Together, 
this gives r c ~ \J~A and r r independent of A. Hence, in 
the static regime, the coherence time, which is then just 
r c , would grow as y/A, but for sufficiently large A, this 
will always take us in the motional narrowing regime, 
t c 3> T r , and the coherence time will be t 2 /t 7 . ~ A. 
This argument shows that the decoherence is a meso- 
scopic phenomenon; as in other ordered phases such as 
ferromagnets and lasers, the finite condensate has a finite 
crgodic time, beyond which the order-parameter correla- 
tions disappear. However, in the thermodynamic limit 
A — > 00, this ergodic time diverges, and true symmetry 
breaking can occur. 

In summary, we have shown that the available ex- 
perimental results are well explained by a model of a 
microcavity polariton condensate as a pumped, dissipa- 
tive system where the main decoherence process is the 
combined effect of number fluctuations and inter-particle 
interactions. However, our model also predicts that a 
regime of motional narrowing should be accessible at 
higher pump powers, which would lead to much longer 
coherence times. This would make the polariton conden- 
sate much more suitable for experiments involving the 
manipulation of its quantum state, such as generating 
Joscphson oscillations, ultimately leading to applications 
in quantum information processing. 
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